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We present detailed numerical simulations of a laser phase stabilization scheme for LISA, where 
both lasers emitting along one arm are locked to each other. Including the standard secondary noises 
and spacecraft motions that approximately mimic LISA's orbit, we verify that very stable laser 
phases can be obtained, and that time delay interferometry can be used to remove the laser phase 
noise from measurements of gravitational wave strains. Most importantly, we show that this locking 
scheme can provide significant simplifications over LISA's baseline design in the implementation of 
time delay interferometry. 
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I. INTRODUCTION 

The LISA mission consists in three spacecraft (S/C) to be launched over the next decade in order to implement 
laser interferometry in space for the observation of gravitational waves of galactic, extra-galactic, and potentially 
cosmological origin. The S/C will be placed on heliocentric orbit in a triangular configuration, where each arm of 
the triangle will be approximately 5 x 10 9 m long. Gravitational waves will be detected by monitoring the relative 
position of the S/C with high accuracy, leading to a gravitational wave strain sensitivity of ~ 10 _23 /\/Hz from 0.1 
mHz to 0.1 Hz. 

One of the most significant technical challenge to achieve this level of accuracy is the cancellation of the laser phase 
noise. A number of signal processing techniques, known as Time Delay Interferometry (TDI), have been developed 
in order to synthesize noise-cancelling interferometers from the various phase measurements that LISA will provide 
[1 H H H IE H While these techniques will certainly be very useful, they impose a number of relatively strict 
requirements on the system 8], such as a good knowledge of the arm lengths at all time, precise synchronization of 
the clocks between S/C, etc. 

An intriguing idea was recently proposed in Ref. @ to significantly reduce the laser phase noise: the very stable 
arms of the LISA constellation could be used as frequency references to lock the lasers. We describe the control loop 
topology proposed by (referred to as "laser self- locking" ) in section^ where we also present a time domain analysis 
of the behavior of the transients introduced by the locking procedure. We then present in section lTTTI thc results from a 
set of realistic numerical simulations of LISA with laser self-locking, where all standard secondary noises are included, 
and where the orbital motion of the S/C is approximately included. We finally show in section Hvl results for the 
use of TDI to bring the reduced laser phase noise below the other noises. We describe in particular how some of the 
technical challenges of implementing TDI in LISA's baseline design are simplified by using laser self-locking. 

II. LOCKING SCHEME 

As illustrated by Figure ^ we use the same locking scheme as although we add realistic acceleration and optical 
noises to our model. Open loop laser noises p®^ and p° 3 L are generated on-board S/C 1 and 3, respectively. They 
then follow a similar path, p®^ is collimated by the telescope on S/C 1, where it acquires noise A31 from the S/C 
motion, before being propagated down arm 2 towards S/C 3. Upon reception after a travel time Zy? jL0| . it is exposed 
to the noise A13 from the motion of S/C 3. It is bounced off the proof mass (acquiring S/C noise 2A13 and proof 
mass acceleration noise 2^13), and interacts with the local laser on a photodiode, where the optical noise 7113 is added 
to produce the electrical signal S13. This signal would be the normal phase measurement in the absence of locking. 
In the present configuration, however, the signal S13 is fed back to the laser phase actuator through a controller G 13 
in order to close the control loop. 
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FIG. 1: Diagram of the control system. S/C 1 is left of the 16.7 delay lines representing arm 2, while S/C 3 is to its right. 



The equations that relate P31 and P13 to the other noises and to the controllers gain are derived in Q . Although the 
authors of [2j note the presence of a strong, slowly decaying quasiperiodic signal with period equal to the round-trip 
light-time, they do not include it explicitly in their expressions for the noise spectra. This fact was first pointed out 
by Tinto [Tl| and others. However, we show below that the quality factor of this signal is very high, so that its power 
spectrum consists only in narrow peaks at harmonics of the inverse of the round-trip light-time. Its effect on data 
analysis should therefore be negligible, just like high quality factor oscillations from wire resonances, for instance, 
have a negligible effect on the data analysis for ground-based interferometers. 

As a simplification of the diagram of Fig. Q neglect all secondary noises, and assume that S/C 3 is an ideal mirror, 
so that 

P3l(t) = P3l(t) ~ GSI * [P3l(t) -P3l(t ~ 2L 2 )} , (1) 

where the * operator denotes convolution. It is assumed that this equation holds for t > 2L 2l with the control loops 
being closed at t = 0. For t < 2L 2 , is is assumed that no light is reflected off S/C 3 back to S/C 1. Without loss of 
generality, the solution to Eq. ^is assumed to take the form 

P3i{t) = f(t)Q(t)+P(t) (2) 

for t > 2L 2 , where 

Q(t) = Q(t-2L 2 ), (3) 

f(2L 2 ) = 1, (4) 

and P(t) is an arbitrary function with no power at frequencies / = n/2L2, for n — 1, 2, ... Substituting into Eq. ^ an d 
rearranging, 

f(t)Q(t) + G 31 * [f(t)Q(t) - f(t - 2L 2 )Q(t - 2L 2 )} 

+P(t) + G31 * [P(t) - P(t - 2L 2 )\ = P 2t{t). (5) 

In order for P to not have any periodic component at frequency 1/2L 2 , the first line in Eq. 0must vanish. Using Eq. 
Eland assuming that the convolution is dominated by the gain g at frequency 1/2L 2 , we have 

[(l+g)f(t)-gf(t-2L 2 )]Q(t) = 0, (6) 

i.e., 

f(t) = -^-f(t-2L 2 ). (7) 
1+3 



3 



The solution to this equation is 



/(*) = e- t/X 



(8) 



where 



A = 



2L 2 



2gL 2 . 



(9) 



log(l + l/g) 



Since g is typically very large, the quasiperiodic transient decays on a large number of round-trip light-times, and 
therefore produces only a few narrow peaks in the noise spectrum, which have a width ~ 1/ A at frequencies that are 
multiples of l/2L 2 . 

It is also straightforward to show that the function Q is mostly defined by the open loop laser noise from both S/C 
in the time L 2 just after the loops are closed. Colloquially, the lasers "load" the arm buffer with ~ 16.7 s of noise 
after the loops are closed. As soon as light reaches the other S/C, the lasers will try to track the phase of the light 
stored in the arm, so that the initial noise effectively becomes a reference signal for future frequency measurements 
in the closed loop configuration. 



According to Fig. ^ we assume that the electronic signal out of the phasemeter on every S/C is filtered by a high 
gain digital controller, and is fed back to the laser using a phase actuator, which frequency response is described by a 
pole at DC (as in This DC pole provides a 1// attenuation, and a phase shift of —n/2 radian at all frequencies. In 
addition, the time delay for the propagation of the light down the arms provides a transfer function with a significant 
attenuation and a +7r/2 phase shift at low frequencies. Its amplitude goes to zero at multiples of the frequency defined 
by the inverse of the round-trip light-time, where the phase also drops to —ir/2. 

In order to have a controller that is stable at all times during the ramping up of its gain, the total phase of the 
phase actuator, time delay, and controller transfer functions must be larger than —n at all frequencies by some phase 
margin <f>. We chose a conservative <j> = 0.35 rad. Well below the first round-trip light-time resonance at 0.03 Hz, 
the phase shift from the delay line cancels the phase shift from the phase actuator, so that n — cj) is available to the 
controller, and a steep roll-off of f~ 8 is possible (/ _1 " 8 with the actuator DC pole). At its first null, the phase of 
the time delay transfer function is —ir/2, so we have to provide +</> of phase compensation at 30 mHz. Given the f~ l 
response from the actuator, the gain then goes as J -0 ' 8 . At higher frequencies, the phase shift resonances are closely 
spaced, so we must also provide +</> of phase compensation at all frequencies and maintain the f~ 8 attenuation. 

This design was used to construct the controller with response shown in Fig. [3 (note that the top two panels in 
Fig. [^include the effect of the DC pole of the phase actuator). The gain was scaled to provide a unity gain frequency 
around 10 kHz, so that this controller provides close to 180 dB of attenuation at 0.1 mHz. Although a bit more gain 
might be obtained by better fitting the optimal phase response (— n + (f) around 10 mHz and below, this controller is 
a fairly aggressive design that allows the investigation of the limits of the self-locking procedure. 



An optimized code was used to simulate the system with realistic parameters. At a sampling frequency of 100 kHz, 
open loop laser noise with an approximate f~ 2 power spectrum was generated by summing successive data from a 
white noise series. At the first time step, each open loop laser noise datum was propagated through the diagram of 
figure^m order to calculate S31 and S13, as if the feedback through the controller was not present. At the second time 
step, the values of S31 and S13 obtained from the first time step, and filtered by the controller response G13 or G31, 
were also added to the open loop laser noise. This procedure was then repeated for all other time steps, effectively 
producing an integration of the system through the equivalent of the Euler scheme. 

The proof mass acceleration noise (Sij) was an approximate f~ 2 noise generated like the open loop laser noise, 
while the optical path noise (ny) and the spacecraft acceleration noise (Ay) approximately had f 2 power spectra, 
and were generated by taking the difference of successive data in a white noise series. At 1 mHz, the ratio of the rms 
values p° L : Ay : <5y : ny was 1 : 4.88 • 10~ 10 : 1.58 • 10" 11 : 3.95 • 10~ 12 . 

A strong anti-aliasing low-pass filter with corner frequency at 1 Hz was applied to the f 2 noises in order to avoid 
the aliasing of high frequency noise back to lower frequencies when the arms lengths were changing. For an arm 
length rate of change L > 0, a given noise perturbation at frequency / is Doppler shifted to a different frequency 



A. Controller design 



III. NUMERICAL SIMULATIONS 
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FIG. 2: Characteristics of the controller. Top: amplitude response of the controller and phase actuator (DC pole) vs. frequency. 
Middle: phase response of the controller and phase actuator. Bottom: phase response of the control loop, including the 
controller, 33.4 s delay line, and the actuator. The two horizontal dashed lines are at ±7r =p 0.35 rad. Circles and crosses 
represent the frequency location of zeros and poles, respectively. 

/' = /(l + L/c) for every arm round-trip. Over hundreds or thousands of round-trips, a high frequency perturbation 
can easily be pushed above the Nyquist frequency, and be aliased back to contaminate the measurement spectrum. 

The propagation delay down the arm was simulated using a data buffer with varying length. At any given time, the 
transmitted value was taken from the linear interpolation between the two nearest points in time. The interpolation 
was critical in order to prevent the introduction of numerical artifacts in the calculations. Ref. [2 shows that the 
motion of the S/C can be modeled by arm lengths that are varying with time in a sinusoidal manner. The peak rates 
of change of two of the three arms are approximately equal (~ 1 m/s), and are significantly less than that of the third 
arm (~ 13 m/s), while the period of the variation of the third arm length is longer that the period of the other two 
arms by a factor of ~ 3 (1 year vs. 4 months). As a rough approximation of these orbits, we assume a constant rate 
of change that is equal to 1 m/s for arm 1, and 10 m/s for arm 2. All arm-lengths are assumed to be 16.7 s at the 
beginning of the simulations. 

A. Results 

Figure shows the amplitude spectra of the laser noise (closed and open loop) for the 10 m/s arm change, together 
with the theoretical expectations for a /~ 2 noise and for the closed loop noise. The closed loop model is constructed 
from 9, Eq. (12)], using the known response of our controller. The fit is excellent, except in the vicinity of the 
resonances at multiples of the inverse of the round-trip light-time (Fig. QJ. At those frequencies, the large quasiperiodic 
transient introduces additional noise. Because the transient is decaying over a rather long time-scale, it is quite narrow 
in frequency, and does not perturb the measurement spectrum much. 

Figure shows as a function of time the envelope of the closed loop laser noise, which is dominated by the 
quasiperiodic transient. The fit to an exponential is reasonably good, with a decay time of 1.49- 10 5 s. For comparison, 
the controller gain at / = 1/(2 x 16.7s) = 0.03Hz is 72.8 dB, so that the predicted transient decay rate from Eq. [5]is 
1.46 • 10 5 s, in good agreement with what is observed in the simulations. 

IV. TIME DELAY INTERFEROMETRY 
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From the diagram of Fig. ^ we can write: 

«3i(*) = P3i(i) + n 3l (t) - {26 31 (t) - A 31 (t) + p 13 (t - L 2 [t])} , 



(10) 
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FIG. 3: (Color online) The open (top lines) and closed (bottom lines) loop laser noise amplitude spectra (arbitrary units). The 
upper straight (red) line is the theoretical model for the f~ 2 noise, while the lower continuous (green) curve is the model for 
the closed loop laser noise. The spectra were constructed from 7 ■ 10 4 s of data sampled at 100 Hz. 
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FIG. 4: (Color online) Close-up of figure[J]fbr 0.02Hz < / < 0.1Hz. 



and similarly for the other signals S13, S12 and s 2 i- Here, S31 and s 13 are the signals measured on S/C 1 and 3, 
respectively, for the phase difference along arm 2, while S21 and S12 are the signals on S/C 1 and 2 for the phase 
difference along arm 3. On a given S/C, the relative phase between the two lasers is measured by bouncing the light 
from one laser off the proof mass on its optical bench, and then sending the reflected light through an optical fiber 
to the other optical bench, where it interferes on a photodiode with the laser from that bench. The measurement 
resulting from sending p%\ to the bench where P21 is located, for instance, is given by 

vii{t) =pai(t) -p2i(t) - 2<5 3 i(*) + 2A 31 (i) + /i 2 i(t), (11) 

where H21 is the noise introduced by the fiber. 

Time delay interferometry (TDI) is a set of combinations of the s,j and measurements that have the property of 
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FIG. 5: The logarithm of the amplitude of the envelope of the closed loop laser noise as a function of time (continuous line), 
and the best-fitted straight line (dashed line). 



cancelling the laser noises (pij) and the S/C acceleration noises (Ay) while preserving a good response to gravitational 
waves. Successive versions of TDI variables of increasing complexity have been developed for LISA configurations 
with static arms 0,0, EJL1- an d for rotating and stretching arms |jE10]- The latter, more complicated versions 
are known as "second generation" TDI variables, and are necessary for the small errors in the cancellation of the 
laser noise to be below secondary noises, when the initial laser noise level is that of a stabilized laser without locking, 
~ 30Hz/VHi. 

Figures HH and [7\ show the first generation TDI variable X that is formed with self-locking. For reference, X is 
constructed in the following manner: 

X(t) = s 12 [t - L 3 {t) - 2L 2 {t)] - a 13 [t - L 2 (t) - 2L 3 (t)} + s 21 [t - 2L 2 (t)] - s 31 [t - 2L 3 (t)} 

+s 13 [t - L 2 (t)} - s l2 [t - La(t)] + sax [t] - sai [*] 

~ {-r 3 i[t - 2L 2 (t) - 2L 3 {t)] + r 31 \t - 2L 2 {t)] + r 31 [t - 2L 3 (t)} - T 31 [t] 

+T 21 [t - 2L 2 (t) - 2L 3 (t)} - T 2X [t - 2L 2 {t)\ - r 21 [t - 2L 3 (t)} + r 21 [t]} , (12) 

where 

L 2 (t) = 16.7s (13) 
c 

and 

L 3 {t) = 16.7s +^^t. (14) 
c 

It should be noted that the time series Sy and used to construct X are only sampled at 100 Hz, and that simple 
linear interpolation between the two nearest points is used to apply the right time-shift. In other words, the fact 
that the data from the simulations and the theoretical model agree very well in Fig. shows that laser self-locking 
allows the construction of X with rather loose requirements on the accuracy of the time-shifts. In fact, we have 
verified directly that adding a fixed delay of 100 /is to L 2 when constructing X does not degrade appreciably the 
noise spectrum. With a 1 ms delay added, however, the X amplitude spectrum is increased by a factor of ~ 10 at 
frequencies above 1 mHz. A fixed delay is used for the whole simulation (a total of ~ 19.4 hours) to simulate the 
case where the arm length measurements are not being performed continuously; this could be the case if, for instance, 
these measurements perturb the gravitational wave measurements, or if they involve the exchange of information with 
ground stations. 
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FIG. 6: (Color online) The amplitude spectrum of the X TDI variable from our simulations (blue), and the theoretical spectrum 
(green) constructed from the proof mass acceleration noises (8ij, all with amplitude spectrum as described by the dotted line 
with a negative slope), and from the optical path noises (ny, all with amplitude spectrum as described by the dotted line with 
a positive slope). The spectrum was constructed from 7 ■ 10 4 s of simulated data sampled at 100 Hz. The amplitude units are 
arbitrary, but are the same as in Fig. [3 
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FIG. 7: (Color online) Close-up of Fig. ©for 0.02Hz < / < 0.1Hz. 



Without laser locking, the arm lengths have to be known to a much greater accuracy in order to implement TDI. 
Ref. H shows that the required accuracy is of the order of 100 ns. We have verified directly by constructing X as 
described above with open loop simulation data that the simulated spectrum was 3-4 orders of magnitude above the 
theoretical expectations. This, of course, is a direct result of our 100 Hz data sampling and two-point interpolation 
scheme, which are far from being accurate enough to implement TDI with a large laser noise. 
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V. CONCLUSION 

Starting from the idea presented in Ref. to lock the two lasers of a LISA arm to each other in order to improve 
their phase stability, we have implemented a realistic simulation of the system, which includes a stable control loop 
design (with unity gain frequency at 10 kHz), significant arm length variations (~ 10 m/s) from the orbital motion of 
the S/C, and the standard secondary noises (optical path noise, proof mass acceleration noise, and S/C acceleration 
noise). We have shown that a reduction of the amplitude of the laser phase noise by as much as seven orders of 
magnitude was indeed possible at low frequencies. We have also shown that the large quasi-periodic transient that is 
produced by the locking procedure is very stable and decays exponentially on a time-scale of hundreds of thousands 
of round-trip light-times. We argue that this is the reason why it only adds noise in narrow frequency bands, so that 
its presence is inconsequential for the data analysis. 

We have also used our simulated data to implement Time Delay Interferometry (TDI). Our results show that the 
same level of noise can be achieved as in the baseline design where no locking is used, except in narrow bands at 
the frequency defined by the inverse of the round-trip light-time, and its harmonics. Since the laser noise is much 
smaller when self-locking is used, the accuracy on the arm lengths required to implement TDI can be significantly less 
stringent than in the baseline design, ~ 100/iS vs. ~ 100 ns. We also verify that the first generation TDI variable X 
can bring the laser phase noise below the secondary noises, even in the presence of significant time dependent changes 
of the arm lengths ("flexing"). 

Our results suggest that laser self-locking could be used to move technical complexity and risk from the implemen- 
tation of second generation TDI to the implementation of the control loops required to lock the lasers. It appears 
that self- locking and the baseline design are mostly equivalent from the point of view of the instrument sensitivity, so 
that the viability of both schemes should be carefully evaluated from an engineering point of view. In particular, it 
should be verified that more detailed orbital models, which might include rotation of the constellation in addition to 
the stretching of the arms, do not lead to degraded performances for the self-locking scheme. 
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